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Abstract 

The principal mechanisms of Auger recombination of nonequilibrium carriers in 
semiconductor heterostructures with quantum wells are investigated. It is shown 
for the first time that there exist three fundamentally different Auger recombi- 
nation mechanisms of (i) thresholdless, (ii) quasi-threshold, and (hi) threshold 
types. The rate of the thresholdless Auger process depends on temperature only 
slightly. The rate of the quasi-threshold Auger process depends on temperature 
exponentially. However, its threshold energy essentially varies with quantum 
well width and is close to zero for narrow quantum wells. It is shown that the 
thresholdless and the quasi-threshold Auger processes dominate in narrow quan- 
tum wells, while the threshold and the quasi-threshold processes prevail in wide 
quantum wells. The limiting case of a three-dimensional (3D) Auger process is 
reached for infinitely wide quantum wells. The critical quantum well width is 
found at which the quasi-threshold and threshold Auger processes merge into a 
single 3D Auger process. Also studied is phonon- assisted Auger recombination 
in quantum wells. It is shown that for narrow quantum wells the act of phonon 
emission becomes resonant, which in turn increases substantially the coefficient 
of phonon-assisted Auger recombination. Conditions are found under which the 
direct Auger process dominates over the phonon-assisted Auger recombination 
at various temperatures and quantum well widths. 

I. INTRODUCTION 

Two recombination processes are predominant in semiconductors at high excitation levels: 
(i) radiative recombination and (ii) nonradiative Auger recombination related to the electron- 
electron interaction. For homogeneous semiconductors, mechanisms of Auger recombination 
(AR) have been extensively studied@~i. In narrow gap semiconductors occur AR processes 
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involving two electrons and a heavy hole (CHCC Auger process) or an electron and two heavy 



these processes are of threshold nature, and the rate of Auger recombination changes with 
temperature exponentially 01. The only exception are semiconductors in which the spin-orbit 
splitting is close to the energy gap (GaSb and InAs). Under certain conditions the rate of the 
CHHS process in these semiconductors depends on temperature only slightlyi. It is commonly 
believed that in weakly doped semiconductors phonon-assisted AR processes are predominant 
at low temperature and high excitation levels!'!. Because of the large momentum transferred to 
a phonon, the threshold for heavy holes is removed and phonon-assisted Auger process becomes 
a power-law function of temperature. 

Single semiconductor heterostructures, quantum wells, quantum wires, quantum dots are 
spatially inhomogeneous owing to the existence of barriers. The presence of a heteroboundary 
affects not only the energy and wave functions of carriers, but also the macroscopic properties 
of heterostructures 0, which is of primary importance. It is commonly believed that the AR 
mechanism in quantum wells is the same as in a homogeneous semiconductor Nev- 
ertheless, the presence of a heteroboundary strongly affects the electron-electron interaction 
in quantum wells, and this influence is of fundamental nature. The heteroboundary lifts re- 
strictions imposed on processes of electron-electron interaction by the energy and momentum 
conservation laws. Namely, the conservation of quasi-momentum perpendicular to the het- 
eroboundary breaks down. In turn, this leads to the appearance in heterostructures of new 
thresholdless channels of Auger recombination!. The rate of the thresholdless AR process is a 
power function of temperature. The first direct experiment aimed at observing the thresholdless 
AR channel at T = 77K was reported inEl. At low temperature the thresholdless process, in 
contrast to that with a threshold, is rather an efficient channel of nonradiative recombination, 
and for narrow quantum wells and high concentrations it dominates the phonon-assisted AR 
process. The existence of a thresholdless matrix element of electron-electron interaction also 
exerts strong influence on the phonon-assisted AR process. The latter becomes resonant and 
is much enhanced as compared with the 3D phonon-assisted Auger process. 

A detailed analysis of the threshold and thresholdless AR mechanisms has been performed 
for a single heterobarrier!. Conditions were studied under which the thresholdless channel 
dominates over the threshold one. For quantum wells no such detailed analysis has been done. 
Taylor et a/H considered the possibility of removing the threshold for AR in quantum wells 
upon transition of excited carriers to the continuos part of spectrum. However, no microscopic 
theory of the thresholdless process was given in this work and no theoretical analysis performed 
of the competition between the threshold, quasi-threshold, and thresholdless AR mechanisms 
at various temperatures and quantum well widths. Only the thresholdless AR channel, cor- 
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responding to small momenta transferred in Coulomb interaction of particles (for the CHCC 
process) with spin-orbit interaction neglected, was considered 

The aim of the present work is to investigate theoretically the principal mechanisms of 
AR for nonequilibrium carriers in semiconductor quantum wells. It will be shown that there 
exist three fundamentally different AR mechanisms: (i) threshold mechanism similar to an 
Auger process in a homogeneous semiconductor, (ii) quasi-threshold mechanism with a thresh- 
old energy strongly depending on the quantum well width, and (hi) thresholdless mechanism 
inoperative in a homogeneous semiconductor. For the threshold AR process in a quantum well 
the threshold energy is close to that in a homogeneous semiconductor. Conversely, owing to the 
small value of the threshold energy, the rate of the quasi-threshold process in narrow quantum 
wells depends on temperature only slightly. For this reason there is no clear distinction between 
mechanisms (ii) and (hi) in sufficiently narrow quantum wells, and they may be considered as a 
single thresholdless AR process. With increasing quantum well width, the threshold energy of 
the quasi-threshold process increases and approaches the bulk value. A totally different behav- 
ior is characteristic of the thresholdless AR mechanism. With increasing quantum well width, 
its rate falls dramatically and, on passing to a homogeneous semiconductor, this mechanism 
ceases to be operative. Conditions will be found under which the thresholdless AR mecha- 
nism dominates the threshold one. The critical quantum well width will be found , at which 
the quasi-threshold and the threshold Auger processes merge and form a single 3D AR process. 
Also, the phonon-assisted AR in quantum wells will be studied. It will be shown that for narrow 
quantum wells the act of phonon emission becomes resonant, which in turn increases substan- 
tially the coefficient of phonon-assisted AR. Conditions will be found under which the direct 
AR process dominates the phonon-assisted AR process at various temperatures and quantum 
well widths. 

II. PRINCIPAL EQUATIONS 

To analyze the AR mechanisms and find the rate of the Auger process, wave functions of 
charge carriers are to be known. As already established for bulk Auger processes, the wave 
functions of carriers must be calculated using the multiband approximation!. We will use 
the four-band Kane's model, the most adequately describing the wave functions and energy 
spectrum of carriers in narrow-gap A///By semiconductors^. 

1. Wave functions in a homogeneous semiconductor 

For most A///By semiconductors, wave functions of electrons and holes in the center of the 
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Brillouin zone are described by the representation for the conduction band and by the Tj 
and Tg" representations for the valence band. Of these the first two and the last are doubly 
and fourfold degenerate, respectively. The corresponding equations for wave functions may be 
written in differential form. Commonly, the basis wave functions of the conduction and valence 
bands are taken in form of eigenfunctions of the angular moment um@'0. However, another 
representation of the basis functions is more appropriate for our purposes: 

\s1),\sl),\xV,\xl),\yV,\yl),\zV,\zl), (1) 

where \s) and \x), \y), \z) are the Bloch functions of s and p type with angular momenta of 
and 1, respectively. The former describe the state of the conduction band and the latter 
the state of the valence band at the T-point. Arrows denote the direction of spin. The wave 
function of carriers ip may be presented in the form: 

V> = *,|s) + *|p), 

where \l/ s and ^ are spinors. In the vicinity of the T-point the equations for \l/ s and ^ envelopes 
written in the spherical approximation follows: 



(E c -E)V e - z&7W = 0, 

(E v -S- E)V - ift 7 W s + £(ti + 47 2 )V(V*)- (2) 
-g (7i - 27 2 ) V x [V x + i8[tr X *] = 0. 

Here 7 is the Kane's matrix element^ having dimension of velocity, 71 and 72 = 73 are the 
generalized Luttinger parameters^, 5 = A so /3, A so is the spin orbit splitting, E v and E c are 
the energies of the lower edge of the conduction band and the upper edge of the valence band, 
m is the free-electron mass, cr = (a x , cr y , a z ) are the Pauli spin matrices. If, instead of using 
the Luttinger parameters, the heavy hole mass describing the interaction with higher bands is 
introduced phenomenologically, then equations (|2|) are transformed into equations derived by 
Suris 0. It is easy to verify that equations (0) are identical to those commonly use d0@-@. In 
the first equation in the system (Q) we neglect the term with mass of free electrons. 

A. Hole states 

The expression for can be found from the first equation of the system (^). Substitution 
of ^ s into the second equation gives: 

h 2 h 2 

-EV + V(W) Vx[Vx*l+i«J[(TX*l = ) (3) 

zrrii 2rrih 

where 
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ml 1 = + m" 1 (7 1 + 47 2 ), ml 1 = m' 1 ^ - 2~f 2 ). 

E g + 5 — E 

Here coincides with the heavy hole mass, and mi with the light hole mass in the case of 
zero constant of spin-orbit interaction; E g = E c — E v is the semiconductor forbidden gap. For 
the sake of convenience, it is assumed that E v = 5. This choice is due to an increase in the 
heavy hole and light hole energies at the T-point by 5 and a decrease in the spin split off (SO) 
hole energy by 25 under the action of spin-orbital interaction (see eq. |7|). Equation (|3|) can be 
simplified by introducing new functions 

= divSP and 77 = crrot^. (4) 

After taking the divergence and rotor of equation ([|), multiplied by a it is transformed into 
a system of two differential equations 

-£0+^ + ^ = 0, 
-(E + 5)77+^^-2^ = 0. 

Fourier transform of these equations gives hole spectra for a homogeneous semiconductor 



E + £—k 2 i5 

2m; 

-2iS E + ^k 2 + 5 



0. (6) 



The characteristic equation has two roots 
E h2 = ---—(mT 1 + m h 1 )± 



^ ^-^(-r 1 --, 1 )! 2 - (7) 



4 v 1 h ' \ \2 4 

It should be noted that mi depends on energy (see eq. (§)). At the T-point (k = 0) we have 
the roots E 1 = 5 and E 2 = —25. The positive solution corresponds to light holes, and that 
with negative sign to SO holes. 

In the vicinity of the T-point the energies Ei t2 can be expanded into a series in terms of wave 
vector to relate the effective light and SO-hole masses m^, m so and the Luttinger parameters: 

E 1 &6 , E 2 ss -25 , (8) 

2m h[ 2m so 

where 

-1 _ 4 7 2 (71 + 272) -1 _ 2 7 2 71 
mu ~ 3E g + m ' mso ~ 3(E g + 35) + m 

An approximate light hole spectrum can be obtained by means of a widely used 4x4 Hamil- 
tonianEi However, the range of its applicability is rather narrow, since commonly mi ~ 0.1 m/, 
and the expansion (^j) is only valid when E -C ^A so . Moreover, such a model cannot describe 
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Auger transitions at all, since the basis states of carriers in different bands are orthogonal. The 
same applies to the spectrum of the SO band. 

The Fourier amplitudes of the wave functions of both light and SO holes can be presented 
in the form (see eq. (|3|)): 

2m h 

where / is an arbitrary spinor related to the previously introduced function by the expression 
<f> = k 2 (E)f. 

The third solution of (|3]) pertaining to heavy holes satisfies the relations div^f = (as a 
consequence \l/ s = 0) and crrot^ = 0. This follows from equation (|5|), since, if = 0, then 
7] = and vice versa. It can be readily seen that 

[<x x * fc ] = -i& h . 

Thus, the dispersion law describing the heavy hole spectrum looks like 

E h = 6-p±. (10) 
The components of the heavy hole wave function must satisfy the equations: 

j ] .7 , <rx * (11) 

k z ^ z t + k x ^ x t + k v fy v f = 

T T v yJ divV = 0. (12) 

k z ^ zi + k x y xl + k y V yl = 

Solving these equations one may obtain the explicit expressions for the wave functions. For a 
quantum well they are given in Appendix A. 



B. Electron states 



In principle, the conventional equations for electrons have the same form as those for holes. 
Since the T-point in the conduction band is only doubly degenerate, and the crystal field 
causes no additional splitting, there is no need to retain terms with parameters Moreover, 
the presence of these terms in the equations for electrons gives a far too exact model. Thus, a 
simplified model will be used for electrons 
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(E c - E)^ s - ih-fV® = 0, 

(E v -8- - ft 7 W s + i5[cr x *] = 0. 



(13) 



The electron energies can be conveniently reckoned from the lower edge of the conduction 
band (E c = 0). This energy will be denoted by £ , so that it would not be confused with the 
full energy of electron E, reckoned from the same level as the hole energy. Introducing into 
equation (|T3|) the functions <fi and r) in the same form as before (see eq. (f|)), one can obtain 



-(E g + 5 + £)<j) + ^A(j) + i5rj = 0, 
-{Eg + 8 + 25)r] - 2i5(p = 0. 

Passing to a Fourier transform, we find the electron dispersion law 

£ £ 2 + £(2E g + 35) + (E g + 35)E g 



k 



E g + £ + 25 



If £ <C E g , S, then the energy is quadratic in wave vector. 



£ 



h 2 k 2 

2m c ' 



where 



27 2 



E g + 26 



(Eg + 35)Eg 

The Fourier amplitude of the wave function is given by 



£ 



h-ik 2 (£) 



id 



£ + E + 25 



[k x (af)} 



(14) 



(15) 



(16) 



(17) 



where / is an arbitrary spinor (see eq. (§)). 



C. Probability flux and the equations near the heteroboundary 



An expression for the probability flux density can be derived from equation (0) by substitut- 



ing; E - -ih& and using then a procedure similar to that employed in quantum mechanics^. 
It can also be derived by the kp method in the second-order perturbation theory. As a result, 
the following expression is obtained in the case of holes for the probability flux density 

= E 9 + 6 '- E m s y* + _ x rot iIr* _ x rot*). (18) 

2ma 2m h 

For electrons in the conduction band this expression takes a simpler form 
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j e = 7[*»** + 



(19) 



The exact procedure for deriving boundary conditions for wave functions at the heteroboundary 
still remains to be devised. However, some approximate methods for solving this problem have 
been developed in recent years. The Kane's parameter 7 usually differs only slightly for AjjjBy 
semiconductors, hence continuity of 7 is usually supposed in literature (see for example0). 
Discrepancy of parameter 7 in a quantum well and barrier region results in small change of 
Auger coefficient (see section IV). Following the method elaborated by BurtEl and asuming 
contunuity of Kane's parameter, we derive from the system (D) Kane's equations which can be 
integrated across the heterobarrier: 

(Eg + 6- E)^ s - zft 7 W = 0, 

-EV - ffi 7 W s + £ V [672 W] + (20) 
l+£A(7i- 27 2 )afc* + ^x*]=0. 

Using these equations and the probability flux density conservation law we derive the boundary 
conditions for the wave-function envelopes. 

2. Carrier states in a quantum well 

The wave functions of carriers in a quantum well may be derived using the symmetry prop- 
erties of the Hamiltonian. Spinless Hamiltonian 7i is invariant with respect to the substitution 
x — > —x. Consider an operator 1Z such that 

K : (x,y,z) -»• (-x,y,z), H = XC VX , (21) 

where I is the inversion operator, and C nx is the operator of rotation by an angle 7r around the 
x axis perpendicular to the plane of the quantum well. 

With account of the spin orbit interaction the Hamiltonian can be written in the form: 

H = H q + ^- 2 [VVx P ]<t, (22) 

where p is the momentum operator and V is the potential energy of an electron in the crystal. 
The last term does not commute with R. Therefore, the symmetry operator T> may be sought 
for as a product of operator TZ and some spin matrix S to be found: T> = 1Z®S. Since inversion 
leaves unchanged the sign of the vector product, the matrix S must satisfy the relations 
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Sa x = a x S 



S 



o~,, 



-o y S , <J X 



"o 


1 " 




"0 -i 




1 


1 







i 




-1 



(23) 



Obviously a Pauli spin matrix a x may be taken for the matrix S: S = a x . 

The functions *&(x, y,z) and P 1 i r (— x, y,z) satisfy the same equation. For this reason the 
eigenfunctions of the Hamiltonian may be sought for as eigenf unctions of the operator T>. 



ty(x,y,z) = uV^(—x,y,z), where v — ±1 



(24) 



The values v — ±1 correspond to carrier states with different symmetry With the wave 
functions chosen in such a way, the boundary conditions may be satisfied at one heteroboundary 
only since at the other they will be fulfilled automatically. Solving equation ( ^4] ) we find the 
necessary conditions for various components of the symmetrized wave function. 



^ sT (x, y, z) = ±$ al (-x, y, z) , ^ xT (x, y, z) = ^ xi (-x, y, z) , 
* v t(s, y, z) = ±V yl (-x, y, z) , ^- A (x, y, z) = ±^ zl (-x, y, z) 
where the sign " + " corresponds to v — 1, and " — " to v — — 1 



(25) 



Corresponding expressions for the components of electron and hole wave functions are given in 
Appendix A. 

III. PROBABILITY OF AUGER RECOMBINATION 

The probability of AR per unit time can be calculated in terms of the first-order perturbation 
theory in electron-electron interaction: 

2tt, 

T 



w, 



\M fi \ 2 d(e f -Ei), 



(26) 



where 



M fi = (*/(n,T" 2 , Vi,V2) 



KOFI - r 2\ 



$(ri,r 2 ) 



(27) 



is the matrix element of electron-electron interaction, r± and r 2 are carrier coordinates, v\ and 
v 2 are spin variables (see eq. (p4|)), e is an electron charge, and kq is the dielectric constant 
of the intrinsic semiconductor, $(ri, r 2 ) is the additional potential arising because of difference 
between quantum well and barrier dielectric constants. The explicit expressions for $(r!,r 2 ) 
are given in the Appendix B. 

Taking into account the antisymmetrized form of the wave functions, the matrix element of 
Auger transition is the following: 
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Mfi = Mi - Mii, 



(28) 



where 

e 2 



Mi = (* 3 (ri,n)*4(r 2 ,i*j) 



+ $(ri,r 2 ; 



*i(n,^2(r 2 ,^)>, (29) 



the expression for M/j may be derived from (|2~9|) by swaping indexes 1 and 2 in the wave 
functions and ^2- Hereafter the indexes I and II in the expressions for the matrix elements 
will be omitted. 

We shall consider two AR processes, CHCC and CHHS, since in fact only these two deter- 
mine the rate of Auger recombination. Strictly speaking such a terminology is inapplicable to 
carriers in a quantum well, since there exists mixing between heavy hole, light hole and SO 
hole subbands. However, as noted above, in the case m c <C the extent of mixing between 
heavy and light holes is low, and the mixing of SO holes with heavy and light holes is negligible 
at A so ^> T. The last condition is nearly always fulfilled for AjuBy semiconductors. For this 
reason we may rely on the above terminology. 



Matrix element of Auger transition 



Evaluations of matrix elements for the CHCC and CHHS Auger processes are similar. For 
the sake of simplicity later in this section the matrix element of the CHCC Auger transition 
will be mainly discussed. However, in the following section approximate expressions for the 
Auger coefficient will be given both for the CHCC and for the CHHS processes. The matrix 
element of an electron-electron Coulomb interaction can be most conveniently calculated using 
a Fourier transform. We take into account that the wave functions of carriers in a quantum 
well are the plane waves along the lateral direction: 

The explicit expressions for wave functions of electrons and holes ipi are given in the Appendix 
A. Then 

Aire 2 1 r°° r a / 2 
M = / / ^(sOV-SCna) x 

Kq Zq J-oo J -a/2 

x (V^ 1 -^ + fax, x 2 , qj) Mxi)Mx2) d Xl dx 2 , (30) 

q = \qx — 54 1 = |g 3 — q 2 \ is the momentum transferred in the plane of the quantum well in 
Coulomb interaction, a is the quantum well width corresponds to the potential $ and the 
expression for it is given in the Appendix B. The integrating over x 2 is limited within the 
quantum well due to the fact that heavy holes, because of their relatively big mass, are usually 
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strongly localized inside the well. Hereafter x denotes the coordinate orthogonal to the quantum 
well plain and p denotes the pair of coordinates in the quantum well plane, q and k are the 
lateral and x- quasimomentum components of particles. 

As it is seen from the equation (|30D the Auger scattering occurs on the one-dimensional 
exponentially decaying potential which depends on the lateral transferred momentum. The 
state of the excited particle may lie in both the continuous and discrete spectrum, with the 
latter situation occurring when the longitudinal momentum of the particle much exceeds the 
transverse momentum. In determining the rate of Auger recombination, both localized and 
delocalized states must be considered as final states of the excited particle. Possibility of an 
electron (hole) transition into a localized or into a free state leads to the existence of different 
AR mechanisms in quantum wells. 

Calculation of the matrix element of Auger recombination for a transition of the excited 

particle into the continuous spectrum. 

For evaluating the matrix element we use the approximation 



where V c and V v are the barrier heights for electrons and holes respectively. Obviously, this 
approximation also implies that k\ + q 2 3> k 2 , i.e. the total momentum of the excited electron 
is much larger than that of the localized one. Integral over x\ coordinate can be found by 
integrating by parts. The n-th antiderivative of the function %l)^e~ qx is: 



V c , V v < E, 



F?(q,x) = (-1) 



n 



(e qx ip 4 (x)Y n ^ 
{k\ + q 2 Y 



e -2qx 



Then the approximate expression for the matrix element Mj may be obtained: 



M » M (1) + M (2) , 



(31) 



where 



Ko(q 2 + kl) 



4ire 2 



( 



F{a/2) jf 



■a/2 



a/2 



e qX2 ijl{x 2 )^ 2 {x2)dx 2 




(32) 



Here 



i 
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F{a/2) = e-^ 2 r is {a/2)^ ls {a/2) 



3V C + V V Kq — K 

v AE g + 



Index s in the ip^ s and ipils implies that only s-components of the wave functions are taken, 
ko is the dielectric constant in the barrier region. 

M = ( 2^,2, / r^t{x)^2{x)^x{x)dx (33) 
K (q Z + fc|) J -a/2 

Note, as the wave functions are spinors the components of ip\ should be multiplied by compo- 
nents of ipx and vice versa the components of ^3 should be multiplied by those of fa- 

In this way it appears that the matrix element of Auger transition splits two parts. The 
first of them is related to the presence of heterobundaries and the second one corresponds to 
the short range Coulomb scattering. The latter can be easily understood as during Auger 
transition a large momentum is transmitted from the localized electron to the excited one 
and this is possible only if the scattering particles find themselves very close to each other. 
Note that both and and, consequently, the matrix element M itself are in fact 

thresholdless matrix elements. Indeed, they are not subject to any restrictions imposed on the 
initial momenta of carriers, k\, k c , kh- However, the mechanisms responsible for the momentum 
nonconservation k\ + k 2 7^ k 3 + k A in the components M^ 1 ' and are different. In M' 1 ' the 
nonconservation is related to carrier scattering at the heteroboundary, and the same mechanism 
gives rise to a thresholdless Auger process in scattering on a single heterobarrier 0. The reason 
why the conservation law breaks down for is that the volume of integration with respect 
to x is restricted to the quantum well region, which results in the appearance of a function of 
the type S * n instead of 5(k). The physical meaning of the above distinctions between the 
matrix elements M^ l > and is that the corresponds to the true thresholdless process 
whose origin is related to momentum scattering by heterobarriers. The matrix element 
corresponds to a quasi-threshold process, and, at the quantum well width a approaching infinity, 
it transforms into the conventional threshold matrix element. 

87re 2 / 3 V c + V V k — k " 



x 



^{q 2 + kl){q 2 + kl) \ AE g k + k 
x (^(fl/2)^(a/2)) (^(fl/2)V 2 (a/2))' (1 ± e"* 8 ). (34) 



Sign ± in the last brackets is chosen according to parity of the production ip^(x)i(j2{x) , + 
corresponds to the even production and — corresponds to the odd one. In the case qa 3> 1 
this exponent may be omitted and the matrix element corresponds to the independent 

scattering at two heterobundaries. Note, that the matrix element M equals to zero if the 
parities of productions ip^(x)ip^(x) and ipX(x)ipi(x) are different. 

Let us now analyze . The integral entering into is proportional to the sum 
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r 4 ^m(x)Mx)Mx) dx^Y. ± sin( t 4 k J a/2 > (35) 

K4 — k 

where k runs through eight different values k = ±k\ ± k2 ± ^3. Of all terms in the sum (|35|) 
the largest is that for which k = k\ + k 2 + k%. i The contributions to the sum from other terms 
are less significant and will be neglected for the sake of simplicity. Then the expression for the 
matrix elements of the quasi-threshold Auger process takes the following form: 

M< 2 >« — c iS h ^ 1 + 2 / 3a A A AA sm(k f -k cl -k c2 -k h )a/2 x 
K (q 2 + kl) Eg 1 + a c f c h ^ _ ^ _ ^ _ ^ 

J q h k c e l<1 ' 2 ' 3 + q c k h , v c = ±v h , 
M . . (36) 

[ q c q h sm0 2 ,3, = =F^h- 

Here 5 is an insignificant phase coefficient, Aj denotes the normalizing constant, v c and 



are the spin indexes introduced according to (|24|), ^2,3 is the angle between lateral momenta 



of the electron and the hole. Note, that the matrix element for transition into continuous 
spectrum (M) was split to and in different way than inHH, in the purpose to make 
the corresponding expressions more clear. 

Calculation of the matrix element of Auger recombination for a transition of an excited 

particle into the discrete spectrum. 

We now turn our attention to analyzing the matrix element of an Auger transition in which 
the high-energy particle remains in bound state ip^. This case corresponds to the approximation 
qA ^> k±. The matrix element can be calculated similarly as earlier: 

/ % 47TP 2 /" a / 2 

M ~ fU^ / WMU>M dx. (37) 

K {q 2 + hi) J -a/2 

This integral can be readily calculated; however, the general formula is rather cumbersome 
and will not be presented here. We shall only make an estimate of M 3 , valid in the case when 
bound carriers are in the ground quantum state. Then we have: 

M (3) « ^—^e isT ^A c A f A c A h a/2aq c q h sin0 2 ,3 {v c = -v h ). (38) 

Here a is a coefficient of order unity, resulting from integrating the product of the envelopes of 
the carrier wave functions over the quantum well: 

f /2 ! A/2/3/4 dar«a/2a, (39) 
Jo 



2 
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where fi = cosfcjX, i-numerates the initial and final states of particles taking part in the AR 
process. For wide quantum wells, in which particles may occupy different bound quantum 
states, a is given by (cf. Q55|): 

Here by the index i is meant summation from 1 to 4, and Oi characterizes the parity of the 
function fi (o"j = 1 and for odd and even functions, respectively). 

IV AUGER RECOMBINATION COEFFICIENT 

To calculate the rate of AR, the probabilities of Auger transition per unit time (^) should 
be summed over all initial and final states of carriers with appropriate weights-occupation 
numbers. 

2n 

G = Y ^ (M 2 )-f l f 2 (l-f 3 )(l-U)S(E 3 + E,-E 1 -E 2 ). (41) 

ki,k2,k 3 ,k4 

Here f\ and f 2 are the occupancies of the initial states and f% and f'4 are those of final states, 

is the squared Auger matrix element, summarized over spins of initial and final states. It is 
more convenient to choose electrons and holes as carriers for the CHCC and CHHS processes, 
respectively. For high-excited states the distribution function / 4 may be set zero. Note that 
instead of 1 — fa, we may write / 3 , where f% is the distribution function for carriers of the 
opposite sign: holes for the CHCC process and electrons for the CHHS process. 

Contributions to the rate of Auger recombination from the matrix elements , and 
can be separated since the excited particles occupy different quantum states. The matrix 
elements and , on the one hand, and , on the other, describe transitions in which 
the excited particle occupies a state of continuous and discrete spectrum respectively. It is 
more difficult to separate the contributions from and M^ 2 \ Even though the physical 

difference between these terms is preserved, there exists a term of interference between them. 
At small quantum well widths the interference is essential, since both processes are in fact 
thresholdless; however, even when the interference is neglected we still obtain a result of correct 
order, reflecting all the main specific features of the AR coefficient as a function of temperature 
and parameters of a structure with a quantum well. For a sufficiently wide quantum well the 
interference between and M^> may be neglected. Indeed, while as a function of 

quasi-momentum shows no extremum, the modulus of exhibits a maximum at the point 
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k&(q) + ^3 = ki + k 2 - When the quantum well width tends to infinity, the maximum at this 
point is of the form of a 5-function. In accordance with the aforesaid, the AR probability for 
wide quantum wells, corresponding to the matrix element , has a maximum (as a function 
of the longitudinal momentum of the heavy hole qh) at higher q^ values than the probability 
associated with With decreasing quantum well width the maxima of these probabilities 

approach each other, and the region of overlapping between these matrix elements becomes 
larger. 

The probabilities of an AR transition within the CHCC process, corresponding to the matrix 
elements AfW and M^ 2 \ are shown in Fig. 1 as functions of the longitudinal momentum of 
the heavy hole at different quantum well widths. It can be seen that the interference between 
the thresholdless process represented by and the quasi-threshold process represented by 
is observed, in accordance with the aforesaid, only for narrow quantum wells. It should 
be noted that the AR probabilities are rather smooth functions of the longitudinal momentum 
of heavy hole (q), since, in calculating them, summation is taken over discrete quantum states 
of carriers. At q close to the maximum value determined by the conservation of longitudinal 
momentum and energy, the AR probability shows a square-root divergence eliminated upon 
integration with respect to q, i.e., in calculating the rate of AR. The probability of Auger 
transition for the CHHS process has a form similar to that for the CHCC process. 

In line with the aforesaid, let us present the rate of AR as follows 

G = G x + G 2 + G 3 , (42) 

where the rate Gi corresponds to a thresholdless Auger process with the matrix element Mi, 
rate G 2 to a quasi-threshold Auger process with the matrix element M 2 , and rate G% to a 
threshold Auger process with the matrix element M3. 

The expressions for the rates G\ and G 2 can be derived from ( f4"T|) by changing in it sum- 
mation over fc 4 by integration and passing from 5-function with respect to energy to ^-function 
with respect to momentum. In what follows we shall study the AR coefficients C related to 
the rate G by 

G = Cn 2 p and G = Cp 2 n 

for the CHCC and CHHS Auger processes respectively. Here n and p are the 2D densities of 
electrons and holes, respectively. For the CHCC process we have: 

^ 327TV ft 7 2 F(A so /E g ) jfrg V c ( W c + V v _ k -R \ 2 / qjkj 1 \ 
X ~ K 2 Ela(a + 2/K c ) 2 (kl + Kl) 2 E g \ AE g k, + £ J \ (gg + k\f h(q h ) / ' 1 ' 

where 
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/ l + 2z/3 \ 2 l + 7x/9 + x 2 /6 . 

r m = — — is a coefficient of order unity, 

K J \ 1+x J (l + x/2)(l + 4x/9) y ' 



Note, that if Kane's parameter 7 is discontinuous than the term -j^F- (J^j should be added to 

( 3 ^ c g Vv ~ kI+rI ) ' wnere Eq c is the electron size quantization energy and £7 = 7 — 7 is the 
difference between Kane's parameters in the quantum well and barrier region. However, this 
addition is usually negligible. The angular brackets in (|43| ) and later denote averaging over the 
heavy hole distribution function. In the case of Boltzmann distribution, which is commonly 
the case for holes, this averaging looks like: 



hn *fc 



J ^ POO hn^h 

(f(Qh,k h )} = ^J2 Qhf(qh,k hn )e " T , where 



Qt 



qr = \j2rrihT /h is the thermal wave vector of heavy holes, kh n is the wave vector corresponding 
to the n-th size quantization level of heavy holes. 

For the CHHS process the following expression for C\ can be derived: 

j,vk fc> c 2 f(^so/e 9 ) n 3 I j4SdMi±^L \ whprP (u) 

1 ~ tcgfi E g ikl + k 2 ) 2 a?{a + 2/k c ) ml {E g - A so ) 3 \ + fcgjs^ + ^) / • ^ 

p.. = (2x + 3(l-x)(l-m so /m h )) 2 l + 2x/3 
2x 2 + (x + 3(1 - x)(l - m so /m h )) 2 1 + x 

In the letter case averaging over distribution functions of two holes occurs. In deriving we 
assumed that E g — A so ^> m h /m so T. 

Similarly, we obtain C 2 for the CHCC process: 

ttV ft 3 7 4 F(A so /E g ) / q 2 c k 2 h + g 2 (£; c 2 + \q 2 c ) 1 - cos (fc / - k h - 2fc c )a \ 

°2 ~ —2 n S IT- I o /.. \ 2 \ /„2 , ;.2\;.. /„ \ L ^7T^ / • l 4C V 



2 

^ 9 



R 5 a(a + 2/ Kc ) 2 \ (g 2 + k 2 h )h(q h ) 2(k f -k h - 2k c f 



Direct calculation of the Auger coefficient C 2 for the CHHS process gives a cumbersome result. 
We present here a simplified expression valid for sufficiently narrow quantum wells at k c ^> q c . 

_ ?r 2 e 4 E c h 3 F(A so /E g ) I 1 - cos(k so - k hl - k h2 - k c )a 

2 ~ AkI E g ml(E g - A) 3 a 2 (a + 2/« c ) \ 2(A; so - k hl - k h2 - k c f 
x Qh2 ((kso + k 2 hl )q 2 hl + q 2 h2 k 2 hl + 2fc 2 1 (q fel q fe2 ) + [q M x q fe2 ] 2 ) \ 
(?M + k 2 i)(q 2 2 + kf l2 )k so I 



And, finally, we have for C3 for the CHCC process: 

32vr 2 e 4 a 1 + \x + ±x 2 1 + § 
/€ 2 ft£ 9 (a + l//t c ) 3 (l + x/3) 2 TT"x \ g 2 (gg, + A; 2 ) 3 



n ~ 32 " 2e4 a 1 + fx + |x 2 1 + §* M g c 2 -f ,\ f47 . 
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Here x = A so /E g , a is a multiplier introduced in (j38|). In the last case we average only 
over discrete quantum states of heavy holes. The threshold momentum q t h is found from the 
conservation law for the energy and the longitudinal component of momentum: 

T—T, _ w , ^(<4 + kl) , h 2 (k 2 cl + k\ 2 ) 



E f (Jkj + qt h )=E g + 



For simlicity we neglected here by longitudinal momenta of electrons, because they are small. 
However, we took into account the size quantization energies of electrons as they change the 
effective band gap in the quantum well. Expanding the energy of excited electron Ef into a 
series in terms of momenta in the vicinity of q t h = Q, where Q is the electron momentum 
corresponding to an energy equal to E g (Q « J im ^.f 9 ) ; we get the following estimation for the 
threshold momentum: 



If the quantum well width tends to infinity, the threshold momentum approaches its bulk value!. 
Account must be taken of the fact that for wide quantum wells with large number of levels the 
introduced multiplier a (see eq. (|5B|)) tends to a 5-function expressing the conservation law for 
the transverse quasi- momentum component: 

or 

a 2 — > T^ a Yj 5 ( k h ± k c i ± k c2 ± kd). 
For large quantum well widths, provided that the condition V c E g holds, the inequality 



C 3 C 2 is valid, since the ratio C 3 /C 2 ~ yV c /E g . Hence, for wide wells C 3 may be neglected 
as compared with C 2 - If V c < E g , the following relation holds for wide quantum wells C 3 /C 2 ~ 



V c /(E g — V c ) > 1. For narrow quantum wells the threshold energy of the CHCC process 
increases (see eq. ([18])) and the AR coefficient (|47D decreases relative to the bulk value by a 



factor 

— £- 3m e E Qc 

g 2l J T g2m h T 

The characteristic width of a quantum well for which this phenomenon becomes essential can 
be readily evaluated from the condition that the exponent is unity: 

2rrih 1 . . 

E c w T- — - ^ a w vr— . 49 
3m c q T 

Thus, at quantum well widths a less than several reciprocal thermal momenta a ^ nqx the 
threshold energy Ef h (a) becomes much higher than the bulk value E^ lh (see Fig. 2). For semi- 
conductor compounds A///By the equality ( f4"9"D is fulfilled at room temperature at a quantum 
well width of order of lOOA. 
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For the threshold CHHS process the heavy hole momenta are not specified by the threshold 
conditions and, therefore, we have to perform integration with respect to them. It seems 
impossible to derive analytically the exact result for the Auger coefficient C3 for the CHHS 
process in view of the fact that the matrix element M3 is rather cumbersome. However, an 
approximate expression can readily be obtained by factoring out the averaged squared matrix 
element from the integrand sign: 

2 2 

2ir 1 r _ q hi+ q h2 

C 3 ~ — (M 3 2 ) / qhidq hl q h 2dq h2 d(f) hl d(f) h 2e "t x 

n Zir^q^ J 

x{ (g A _ esHi+s^ + ^ + (50) 

V 2m so 2m h 2m h J 

Here E g = E g + E 0c + 2E h — E 0so , where E 0c , E Q h E 0so are the energies of size quantization for 
electrons, holes, and SO holes, respectively. Let us introduce a threshold momentum, setting 
it equal to 

2 _ 2(E g - A)m so 

h (2 - fi so ) 

where \i ao = m so /m h . Then the expression for C 3 takes the form 

C^^e^m) (51) 
Taking into account that Q t h 3> k so , we obtain: 



, ^ 256m, o7 rV Q\ h V c kj (1 - A, Q ) 2 
^ h\l {Q 2 th + 2klYE g kl + ^ c 1 + 2X1 ' ' 



where a is a multiplier defined similarly as for CHCC process (see eq. (|38|)), A so is derived from 
A; (see eq. (|A^)) by substitution k so instead of k\. Let us consider in more detail the Auger 
recombination coefficient C2 for the quasi-threshold CHCC process. At a — > 00 a substitution 
can be made in the function being averaged in (|45|): 



l-™(k f -k h -2k e )a _ na _ 

2(k f -k h -2k c ) 2 2 v ' ' K ' 

This formula clearly shows the occurrence of a threshold in this limit, and the coefficient C2 
transforms into a 3D expression on being multiplied by a 2 . For comparison we present both 
the result of Gelmonti for C^d and our limiting expression. 



C ™ = 6 ^ K 2 F 5/2 T1/2 1/2 3/2 6 ^ T > (W) 



16^27^ e^rrich 3 1 _^Em. 

"27 «§ eYt^3\4 



C 2 • a 2 = 6— -i ^ 5/2 _ 1/2 3/2 e »* - . (55) 
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The factor 4 in fl55|) results from the necessity to take into account, in calculating M 2 in 
accordance with (j36|), not only the terms with k = k c i + fc C2 + k^, but also those with k = 
k c i — k C2 + kh, k = —k c i + k C2 + kh, and k = —k cl — k C2 + k h . When the quantum well 
width tends to infinity, all the four terms make equal contributions to C 2 . As can be seen, 
the only difference between the expressions fl54|) and (|55|) is in a numerical coefficient. A small 



discrepancy of ~ 2/3 times is due to the necessity to distinguish between the size quantization 
momenta of the two localized electrons: k c i ^ k c2 , besides the expression (Q) was derived 
when A so tends to infinity, while obtaining the expression (|55|) it was assumed that A so < E g . 
In calculating (|55| ) we neglected the quantity V c as compared with E g . In the general case, 
{C 2 + C3) a 2 should be written instead of C 2 a 2 to make expression (|55|) valid. However, the 
limiting transition from the quasi-threshold to the threshold Auger process (see eq. (^)) can be 
realized only for very large quantum wells. Analyzing the probability of Auger transition as a 
function of the heavy hole momentum, one can obtain a qualitative criterion for this transition. 
As mentioned above, the quasi-threshold Auger process dominates in wide quantum wells. The 
probability of this process has two characteristic extrema (see eq. (f45|)). The first of them 
corresponds to the maximum of squared Auger transition matrix element in the vicinity of the 
threshold momentum. The width of this maximum is of the order of the inverse quantum well 
width. The second extremum lies in the vicinity of the thermal momentum of heavy holes. The 
Auger coefficient C 2 can be estimated as a sum of the probabilities taken in these extrema, 
multiplied by the corresponding widths. Then 

C 2 « C\\Q h « q th ) + C?(Q h w q T ), (56) 

where Qh is the value of the heavy hole momentum: Q\ = k\ + q\\ 

^~ Ae » /T ^. (57) 



Cf a \E th , 

Here Xe 3 ~ 2n/q t h is a characteristic wave length of an electron having energy close to E g . 
Comparing the terms C l 2 h and C 2 one can obtain the criterion of the transition from the quasi- 
threshold to the threshold Auger process: 

/ T \ 3 / 2 E th 

a 3> a c , wherea c = Xe g ( — — ) e T . (58) 



tii 



For semiconductors with an energy gap of the order of 1 eV the critical width (a c ) may be as 
large as several thousand Angstroms. However, the value a c is considerably larger than the free 
path length of carriers in semiconductors. This obviously shows that the correct derivation of 
Auger rate in homogeneous semiconductors should involve momentum scattering process if the 
critical width a c there exceeds the free path length!!. 
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With decreasing quantum well width, the maximum of the probability W2 as a function of 
the heavy hole momentum shifts to lower values (see Fig. 1). This reduces the threshold energy 
of the process, and, as a consequence, makes weaker the temperature dependence of the AR 
coefficient. 

Figure 2 shows the threshold energy of the CHCC process as a function of quantum well 
width for all the three mechanisms of Auger recombination Ci, and C3 separately and for 
the overall process C = C\ + C 2 + C3, found from the formula: 

^(T)=T 2 ^, . = 1,2,3 (59) 

The threshold energy for the quasi-threshold Auger process (see Fig. 2) is less than its bulk 

o 

value. The reason is that the critical quantum well width a c ~ 2000 A is greater than the 
maximum width shown in the figure. The threshold energy for the thresholdless Auger process 
decreases with quantum well width and becomes negative at a certain thickness. This is due to 
the fact that the Auger coefficient C\ decreases with increasing temperature for wide enough 
quantum wells (see Fig. 5). With increasing quantum well width, the total threshold energy 
tends to its limiting value Efy denoted in the figure. 

We now turn our attention to the thresholdless Auger process. As already noted, the prob- 
ability of a thresholdless Auger transition shows no extrema as a function of the heavy hole 
momentum. Therefore, the coefficient C\ has a weak nonexponential temperature dependence. 
This phenomenon was first studied by Zegrya and Kharchenkoi. In addition, the function 
C\{T) is nonmonotonic and has a maximum. The presence of this maximum can readily be 
explained. At low temperature and, correspondingly, small longitudinal momenta of carriers 
their wave functions are nearly orthogonal and the C\ value is small. With increasing tempera- 
ture, the characteristic momentum transferred in Coulomb interaction (approximately equal to 
the thermal momentum of heavy hole) grows. This is the reason why at low temperature the 
Auger coefficient increases with temperature. As the temperature is elevated further, the AR 
coefficient G\ passes through a maximum and starts to decrease, since the long-range Coulomb 
interaction responsible for the Auger process is low for large transferred momenta. The tem- 
perature at which the maximum occurs can readily be evaluated by equating the energy of size 

t; 2 2 

quantization of holes to temperature: T « 2n Za 2 • Note that there would be no such maximum 
if the overlap integral I c h were taken to be proportional to the momentum transferred. Such an 
assumption, having in our opinion no justification for the majority of structures investigated, is 
frequently used (see, e.g., ref.[H]) and gives incorrect expressions for the rate of AR and incorrect 
dependence of this quantity on temperature and quantum well parameters. 

The AR coefficient C\ depends rather strongly on the quantum well width a. Depending 
on which term predominates in (fj3p, the coefficient C\ decreases with increasing a either as 
1/a 3 or as 1/a 5 or as 1/a 7 . In any case, even being multiplied by a 2 , C\ remains a decreasing 
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function of quantum well width. For this reason such a process may be dominant only for 
sufficiently narrow quantum wells. At a pa 1/k c the coefficient C\ exhibits a maximum related 
to the weak overlapping of carrier wave functions. With the quantum well width decreasing 
further, the rate of the thresholdless Auger process falls gradually. The similar expression for 
C\ in the CHCC process was obtained by Dyakonov and Kachorovskii, and Zegrya et al. 

Figure 3 shows the AR coefficients C\ and C2 for the CHCC process as functions of quantum 
well width at different temperatures for a model structure based on InGaAsP. It can be seen 
that, firstly, both C\ and C2 show a sharply pronounced maximum. Interestingly, the positions 
of these maxima are practically temperature-independent. Secondly, the relative role of the 
quasi-threshold process becomes more important with increasing temperature. The threshold 
process is not shown in the figure, since at the quantum well widths considered the coefficient 
C3 is much smaller (by several orders of magnitude) than C\ and Ci- For this reason the 
dependence of C3 on quantum well width is shown separately in Fig. 4. Note that for this 
process the maximum is achieved at a much wider quantum well than for the quasi-threshold 
or thresholdless process. This is in the first place due to the reduction of the threshold energy 
of the threshold process with increasing quantum well width (Fig. 2), rather than to the 
overlapping of the wave functions. The coefficients C\ and C2 as functions of quantum well 
width for the CHHS process qualitatively agree with the curves presented in Fig. 3. 

Figure 5 shows the temperature dependence of the overall AR coefficient for the CHCC 
process and the partial contributions from the thresholdless and quasi-threshold mechanisms 
at different quantum well widths. It can be seen that at low temperature and sufficiently wide 
quantum wells the thresholdless Auger process predominates (C*i > C2), and at high temper- 
ature, conversely, the quasi-threshold process becomes more important (C2 > C\). Therefore, 
the curve describing the temperature dependence of the overall AR coefficient has a character- 
istic shape with a maximum and a minimum. With increasing quantum well width, both the 
maximum and the minimum of the AR coefficient shift to lower temperature and, in the limit of 
infinitely wide quantum well, disappear. Thus, in the case of a homogeneous semiconductor the 
AR coefficient is a monotonic function of temperature. Note that the Boltzmann distribution 
of carriers was used in calculating the Auger coefficients as functions of temperature. At low 
temperature both electrons and holes rule described by Fermi-Dirac distribution func- 

tion. Thus, the average momenta of particles participating in the Auger transition depend on 
temperature only slightly. As a result, at low temperatures the Auger coefficient is a smoother 
function of temperature than in the case of the Boltzmann statistic and it does not tend to zero 
at T — > 0. Figure 6 shows the thresholdless Auger coefficient C\ versus temperature at various 
Fermi energies for quantum wells with different widths. Essential discrepancies between the 
values of the Auger coefficients for the Fermi-Dirac and Boltzmann distributions take place only 
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in the case T -C Ep, where Ep is the Fermi energy for holes. This condition can be realized 
only at very low temperature where the Auger process ceases to be an important mechanism 
of recombination. 



V PHONON-ASSISTED AUGER RECOMBINATION IN QUANTUM WELLS 

At high temperature the threshold AR process predominates in the homogeneous semiconductor 
(C3D oc e t j. However, at sufficiently low temperature such a process becomes exponentially 
weak. In this case the rate of the Auger process is no longer determined by the scattering of 
two free electrons. Mechanisms leading to threshold removal are to be taken into account. It 
is commonly believed that the primary mechanism of this kind is emission or absorption of a 
virtual optical phonon. At the expense of a large momentum transferred to the phonon, the 
AR threshold for heavy holes is removed, and the rate of such an Auger process is a power-law 
function of temperaturei'i'izl. The probability of the phonon-assisted AR is calculated in terms 
of the second-order perturbation theory in electron-electron and electron- phonon interaction^. 
However, the possibility of removing the threshold through interaction with phonons is not 
the only one. At high hole concentrations, hole-hole scattering may become a more effective 
mechanism of removing the threshold for the Auger process. This is why the question of the AR 
mechanism in homogeneous semiconductors at low temperature still remains open. We shall 
discuss this problem in more detail elsewhere§!. In the present work we follow the commonly 
accepted viewpoint that there exists a competition between the phonon-assisted and direct AR 
processes. For quantum wells the situation differs strongly from the 3D case, owing to the 
presence of a thresholdless process even in the first order of the perturbation theory. Therefore, 
it is a priori evident that the conditions under which the phonon-assisted AR dominates over 
the direct process depend strongly on quantum well width. 

As already noted, in the 2D case there exist three AR processes: threshold, quasi-threshold, 
and thresholdless. It will be recalled that in sufficiently narrow quantum wells thresholdless 
and quasi-threshold AR processes are predominant at all temperatures. Below, it will be shown 
that the electron-phonon interaction may affect significantly the rate of AR at low temperature. 
The AR coefficient for a phonon-assisted CHCC process with a threshold matrix element of 
electron-electron interaction is comparatively easily calculated wheni 

E g ^>2fiE g ^>hcu LO ,T (60) 

where uolq is the optical phonon frequency, fi = m c /mh. It can be shown that the coefficient 
of phonon-assisted AR is related to the previously calculated AR coefficient (|55f) by 
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where k = , ftoo is the high-frequency permittivity of the medium, g(a,kth) is a factor 



reflecting the 2D character of holes 2 - 8 ' 

1 1 
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For comparison, we present an AR coefficient calculated for the phonon-assisted 3D case 
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It can be seen that the results for 2D and 3D phonon-assisted AR processes with threshold 
matrix elements of electron-electron interaction are closely allied. A significant difference for 
the case of narrow quantum wells is that the threshold energy E 2 ® increases owing to the 
presence of carrier size quantization levels (see eq. (0)). Correspondingly, the criterion for 
predominance of the phonon-assisted AR process (C^ h ) over the direct threshold Auger process 
(C3) in quantum wells is met at somewhat higher temperature than in the 3D case. However, 
as already noted (see Section III), the rate of the threshold Auger process in narrow quantum 
wells is in itself by several orders of magnitude lower than those of the thresholdless and quasi- 
threshold processes (C3 {C\, C%j). Hence, the phonon-assisted AR process with a threshold 
matrix element of electron-electron interaction cannot compete with thresholdless and quasi- 
threshold processes either (C^ C%, C 2 ). 

Let us now consider a phonon-assisted Auger process with the thresholdless matrix element 
(M ee ) for the CHCC process. Direct calculation of the probability of this process results in a 
singularity at the energy of a hole in the virtual state energy equal to the sum of the energies 
of the optical phonon and the hole in the final state: 
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(64) 



where E s is the energy of the virtual hole, andM ep is the matrix element of scattering of the 
virtual hole by an optical phonon, with the signs plus and minus corresponding to phonon 
emission and absorption, respectively. Note that in the 3D case the summation over interme- 
diate states is in fact reduced to taking a single summand. In the 2D case, nothing of the kind 
occurs, and the sum is taken over discrete quantum states. To eliminate the divergence in the 
denominator of the expression (0), account must be taken of transitions not into stationary, 
but quasi-stationary states, i.e., states with complex energy. In this case the pole (|6"4]) will 
transform into a region of complex energy values: 
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where T = H/t. The characteristic lifetimes r corresponding to these states may vary over a 
wide range, depending on temperature, free carrier density, etc. It only makes sense to consider 
a resonant phonon-assisted process in terms of the second-order perturbation theory when the 
halfwidths of the quasi-stationary hole and the phonon states are less than the energy of the 
optical phonon (Hujlo)- Otherwise, the Auger coefficient must be calculated in the first order 
of the perturbation theory, using the Lorentz function (f(AE) = ^ AE 2 +r2 ) instead of the 
5-function expressing the energy conservation law. For a phonon-assisted AR process with a 
quasi-threshold matrix element of electron-electron interaction both the resonant and virtual 
Auger processes are possible, with the former predominant in narrow quantum wells and the 
latter in sufficiently wide quantum wells. 

In the general case the Auger coefficient for a phonon-assisted process with the quasi- 
threshold and the thresholdless matrix elements M ee may be written as 



C ph = CL + C 2 (65) 



where 
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1 

X (1 + 5m,n) ((™ + ^) 2 7T 2 + Q 2 « 2 ) + (m - n) 2 7T 2 + Q 2 d 2 ' 

The function J n , m (Q) has been calculated for a nondegenerate band. In the case of phonon 
scattering by heavy holes its value will be somewhat lower. However, this fact is insignificant 
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for our purposes. For momenta of bound electrons in the matrix element of electron-electron 
interaction in fl67j ) should be substituted their mean thermal values. The Auger coefficients 
and Cp h correspond, respectively, to phonon emission and absorption. Irrespective of the type of 
the matrix element of Coulomb interaction, the phonon-assisted Auger process is thresholdless. 
This corresponds to the main contribution to the Auger coefficient C p h coming from the hole 
momenta of the same order of magnitude as the thermal momentum. Therefore, in calculating 
C p h me may substitute for the longitudinal hole momentum its mean thermal value. This 
simplifies considerably the expression for the Auger coefficient: 

,2 e Ku LO /T d 2Q I /Q Ea \-\ 



x 7Tcoe^ e ,UJJLUI± ^ r / / at, 



M ee (n, q h + Q| 2 p m 

2 Jn,m{Q)^rzr> ( 67 ) 



2a i m z h 2m h 2m h J 
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|M ee (n, q h + Q| p m 

2 JnMV)"! ( 68 J 



where 



p m = e ? t , g h 



2 V ft 2 

Expressions flS7|) and (B3) can be analyzed easily when the temperature is much lower than 
the optical phonon energy. In this case the thermal momenta of holes q^ may be neglected 
as compared with the phonon momentum Q which is approximately equal to the momentum 
of virtual hole. It is readily seen that the probability of Auger transition as a function of 
Q for phonon emission has two extrema. The first of these corresponds to the minimum 
of the denominator in ( p7D and is related to a resonant Auger transition. Note that for an 
Auger transition with phonon absorption no extremum of this kind is observed and there is no 
resonant process. The second extremum corresponds to the maximum of the squared matrix 
element and, as a rule, is related to a virtual Auger transition. For sufficiently wide quantum 
wells the matrix element of the electron-electron interaction as a function of the heavy hole 
momentum has a form close to the 5-function. In this case the probability of the transition is 
the highest near the threshold momentum, and the process of scattering by phonons is virtual. 
With decreasing quantum well width, the 5-function broadens for the quasi-threshold matrix 
element, and, in addition, the role of the thresholdless matrix element, only slightly depending 
on Q, becomes more significant. This enhances the resonant Auger transition and weakens the 
virtual process. For narrow quantum wells the matrix element of Coulomb electron-electron 
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interaction depends on Q only slightly, and, therefore, the resonant process is predominant. It 
can readily be shown that in this case the following estimation is valid for the AR coefficient 
of the phonon-assisted Auger process: 

Ld LO e 2 m h a 2ir3k(E g ) 2 

Cph * -m^ J ^ m T^T |Mee(go)l ' (69) 

where Q = ^J 2mh £ LO , Hence immediately follows that the phonon-assisted to direct AR 
coefficient ratio has the form 

C p h T hp (M ee (Q )) 2 

(M te Mf • (70) 

where C = C\ + C 2 is the Auger coefficient for the direct process, Y hv — — , r hv is the time of a 

^ T hp 

hole scattering by an optical phonon; is the thermal momentum of holes. It can be seen that 
the phonon-assisted Auger process may dominate over the direct one only in the case when the 
the values of T p h and T are close to each other or at extremely low temperature when the ratio of 
the matrix elements taken at the momenta Qo and qx is large. Note that at high nonequilibrium 
carrier densities, when Auger recombination becomes at all significant, the hole- hole scattering 
turns out to be, as a rule, a much more effective mechanism of relaxation than the hole-phonon 
scattering. This results in a small ratio T ph /T. Therefore, the direct Auger recombination 
dominates over the phonon-assisted process down to very low temperatures. Figure 7 shows 
the coefficients of a phonon-assisted Auger transition (C phl and (C p h2)) corresponding to the 
thresholdless and quasi-threshold matrix elements as a function of temperature for different 
quantum well widths. As a halfwidth G is taken a characteristic value of 20 meV. 



VI DISCUSSION 



Our analysis has shown that for the CHCC and CHHS processes in semiconductor struc- 
tures with quantum wells there exist three AR mechanisms: threshold, quasi-threshold, and 
thresholdless. The coefficients of these processes C±, C2, and C3 depend in different ways both 
on temperature and quantum well parameters: the heights of heterobarriers for electrons and 
holes (V c and V v ) and the well width (see Figs. 3-5). In the limit a — ► 00 the sum of the 
quasi-threshold and the threshold AR coefficients multiplied by squared quantum well width, 
C2 a 2 + C3 a 2 , approaches the bulk AR coefficient C 3D , with the product C\ a 2 approaching zero 
(Fig. 8). For sufficiently narrow quantum wells the 2D AR coefficient multiplied by a 2 exceeds 
the 3D value, owing to the predominance of the thresholdless and quasi-threshold AR processes. 
Thus, Auger recombination in quantum wells is enhanced as compared with that in a homoge- 
neous semiconductor. This enhancement is more pronounced at low temperature. Under these 
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conditions the 3D AR coefficient C 3D is small because of the presence of an exponentially small 
multiplier (see eq. (p4[)). Note that the entire analysis of the AR coefficients (Ci, C 2 , C 3 ) as 
functions of temperature and quantum well parameters is qualitatively applicable to the same 
extent to both the CHCC and CHHS Auger processes. However, since no model structures 
with quantum wells have been specified, we illustrated these relations by the example of the 
CHCC process. 

Note that the AR in quantum wells may be suppressed substantially if the following condi- 
tions are met (V c , V v ) > E g and E 2 — E 1 > E g (Ex and E 2 are the energies of the first and second 
carrier size quantization levels ) il; i.e., in the case when the energy of an excited particle is 
insufficient for a transition into the continuous spectrum or to a next size quantization level. 
For these conditions to be fulfilled, a structure is to be created with deep and narrow quantum 
wells for both electrons and holes. The presently existing technologies can produce structures 
of this kind based on InAs/AlSb § and InAs/GaSb/AlSbi. In these deep quantum wells only 
the threshold AR mechanism, corresponding to the coefficient C3, is operative. This coefficient 
may be smaller by several orders of magnitude than the Auger coefficients for the thresholdless 
and quasi-threshold processes (Ci, C 2 ) in shallow quantum wells ((V c , V v ) < E g ). 

It is also shown that the phonon-assisted AR process undergoes significant changes. Simi- 
larly to the direct AR there exist three different phonon mechanisms (Cp h , Cp h , Cp h ) correspond- 
ing to the threshold, quasi-threshold, and thresholdless matrix elements of electron-electron 
interaction. The first process is quite similar to its 3D counterpart. However, for narrow quan- 
tum wells this process is much weaker than the thresholdless and threshold Auger processes. It 
is this process with the participation of phonons that is considered in the literature to be prin- 
cipal AR process in quantum wellsHS. AR phonon-assisted processes with the quasi-threshold 
and the thresholdless matrix elements of electron-electron interaction may be resonant. At 
low temperature they can compete with direct AR processes. However, owing to the lack of 
any strong temperature dependence in the latter, such a competition is possible at much lower 
temperatures than in the 3D case (see Fig. 7). Thus in narrow quantum wells the direct 
(thresholdless) Auger recombination dominates over the phonon-assisted process in a wider in- 
terval of temperatures than in the 3D case. With increasing quantum well width, the resonant 
scattering by phonons becomes weaker, and we pass to the conventional 3D conditions. 

It should be emphasized once again that at high densities of nonequilibrium carriers in a 
homogeneous semiconductor the phonon-assisted AR process may be less intensive than the 
Auger recombination, with the subsequent hole-hole scattering eliminating the threshold H. 
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VII CONCLUSIONS 



The following principal results were obtained in the work 

1. It is shown that three main mechanisms of Auger recombination are operative simulta- 
neously in quantum wells: 

(a) threshold 

(b) quasi-threshold 

(c) thresholdless 

The first two pass, in the limit of infinitely wide quantum well, into 3D Auger recombi- 
nation, and the coefficient of the thresholdless Auger process tends to zero. 

2. It is demonstrated that in the case of narrow quantum wells the Auger coefficients of the 
quasi-threshold and thresholdless processes show weak power-law dependence on temper- 
ature. Their values much exceed the 3D coefficient related to the squared quantum well 
width. At the same time the coefficient of 2D threshold Auger recombination (a) has a 
higher threshold energy than in the 3D case > Ef®), with the corresponding AR 
coefficient smaller than the 3D expression divided by the squared quantum well width. 

3. It is shown that the quasi-threshold Auger process predominates in wide enough quantum 
wells. The threshold energy of this process is an increasing function of temperature (see 
Fig. 2). The critical quantum well width at which the threshold energy for the quasi- 
threshold Auger process becomes equal to the bulk value, is found. The critical width a c 
strongly (exponentially) depends on temperature. At room temperatures its value may 
be up to several thousand Angstroms for semiconductors with E g « 1 eV. 

4. In narrow quantum wells the phonon-assisted AR process is resonant. At high carrier 
densities it is less intensive than the direct thresholdless Auger process down to very low 
temperature. 
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APPENDIX A: WAVE FUNCTIONS OF CARRIES IN A RECTANGULAR 

QUANTUM WELL 



1. Holes 



Selecting the coordinate system so that the longitudinal component of the wave vector 
coincide with the y axis and performing a Fourier transform in this plane we obtain the following 
expressions for the wave functions of carriers. 

Heavy holes: 

/ 



*l>h(q,x) = Hi 



q cos khX £ 



— ikh sin khX £ 
i —k h sin khX £ + q cos k-^x r\ J 



+ 



(Al) 



+ H 2 



q sin khX rj 
ikh cos khX i] 
i — q sin khX £ — kh cos khX r\ J 



Where f = l/y/2 ( 
Light holes: 



l/y/2 1 } ), H x and# 2 are the normalizing constants. 



i/>i(q,x) = L\ 



+ L 2 



^ ki sin k\x rj — Xiq cos k[X £ \ 

—iq cos k\x rj + iXiki sin k\x £ 
\ —Xih sin k[X £ + Xiq cos kix rj y 

^ — ki cos kix £ — Xiq sin kix r] \ 

—iXiki cos A;/X rj — iq sin fc^x £ 
i — Xiq sin fc^x £ — Xiki cos /c^x 77 y 



(A2) 



si 



ihjjkf + q 2 ) 
E + 6-E 



\L\ cos rj + L 2 sin fc/x £]. 



(A3) 



A, 



E + 25 — h 2 k?/2m h '' 



The wave functions of SO holes are similar to those of light holes. 

A change to functions of another symmetry in the above expressions can be performed by the 
formal substitution £ <-> r\ for \s)-, \x)-, and \y)- components and £ <-> — for |z)-components. 



The wave functions of carriers in the barrier region may be obatined similarly to ( |Al| - |A3|) . 
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If wave functions of two or more particles are considered together, then in general their 
^-components of quasi-momentum cannot become zero simultaneously. A change to a function 
with arbitrary direction of quasi-momentum can be performed using the rotation matrix: 



D^ = R w (g) S t 



(A4) 



where R^ acts on the coordinate components of the wave function, and on the spinor 
components. The Euler angles for a rotation in the plane yz by an angle tp are 



$ = -tt/2 ,Q = (p = 7r/2. 



Thus: 



1 







1 



cos tpj '2 —i sin tpj 2 
i simp/ 2 cos tp/ '2 



(A5) 






costp sin<^ 
— sin ip cos ip 

If the vector q has components q (0, cos tp, sin tp) in the coordinate system x, y, z, then the wave 
function can be written in the form: 

^ q = W = D-^ipo ■ (A6) 

The previously found wave function has a zero subscript. The wave function of heavy holes, 
found using (|A6|) , is written below, as it will be used later. 

q cos khxe~ 1 ^ £ 
-ikh sin khX £ — q cos khX sin <fi r\ 
-kh sin k h x £ + q cos k h x cos <fi rj 

(A7) 



ip h (q,x,<p) = Hi 



+ 



q sin khxe 1 ^ r\ 
+ iJ 2 cos khX rj + g sin k^x sin £ 

— k h cos A;/jX 77 — q sin A^rr cos £ 

The boundary conditions for hole wave functions can be derived by integrating Kane's 
equations across the heteroboundary (see Section II, paragraph l.C). We consider generalized 
Luttinger parameters 71 and 72 also to be continuous for the sake of simplicity. Taking into 

^> m^ 1 we obtain continuity conditions at the heteroboundary for 



account that m l 1 



27 2 



E„+S-E 



the following quantities: 



2 ^ 



(A8) 



3) 



E + 5-E 



divip . 
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Generally speaking the wave functions of holes in a quantum well are superposition of three 
subbands of the valence band: of heavy, light, and SO holes. However, the last subband 
strongly, exponentially, decays away from the heteroboundary with a exponent k so ~ J |. 
As a consequence, this branch mainly affects the derivative of the wave function near the heter- 
oboundary, and its influence on the wave function itself is negligible. It should be emphasized 
that such an approximation is not equivalent to using the 4x4 Hamiltonian from the very 
beginning. We shall seek the wave function as a superposition of the heavy- and light-hole 



~ mi 



This 



subbands. Near the upper edge of the valence band the parameter |A so | 
means that only the first and the third of the boundary conditions ( |A8| ) are applicable. In this 
approximation, light and heavy holes do not mix with each other and have different spectra. 
The heavy hole spectrum coincides with the quantum mechanical spectrum of a particle in a 
rectangular quantum well. For states with even and odd \x > -components of the wave function 
of heavy holes the dispersion equations are following: 



for even states 



tan khCi/2 = — - 
h 

kh 

cot khCi/2 = — for odd states 



(A9) 



For light holes the states with different parities cannot be separated, and the dispersion equation 
becomes more cumbersome: 

,2 



E g + 5 + V c 



Ekf 



x 



E g + 5-E 
E g + 5 + V c 



Kf-q 



T , i , 2Aj - 1 
— k,[ cot kia/2 + k[— 



2A,-1 



x 



Eg + S 

2A,-1 



E 



E kf + q 2 2A, 
-Ki tan kia/2 — ki^~- 



Q- 



2\i - 1 



(A10) 



+ 



E g + 5 + V c - E kf + q 2 



E a + 8-E 



(I 1 



_2X l - 1 

Here Ki and Kh denote the moduli of x quasi-momentum components of light and heavy holes 
in the barrier region respectively, 



A; 



A, 



U v + E + 25 + h Z Kf/2m h 



Note that at q = the light hole states are also split into states with different parities. The 
constants Hi and Lj in (|A1|), (|A2|) are determined by normalization conditions. In particular: 



Hi 



1 



1 



1 2 + k 2 h Ja + 



*h <r+k% 



The opposite is the case for SO holes. The components of the wave functions of light and 
heavy holes oscillate rapidly, and the contribution from them to overlap integrals is negligibly 
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small. Similarly, it is easy to verify that for the so-component ip x and divip/(E g + 8 — E) are 
to be considered continuous. The type of wave functions of SO holes is similar to that for the 
light holes (|A2p,(|X3D. Strictly speaking, with the condition E g — A > U v fulfilled the spectrum 
of spin-split holes is continuous. However, when the rapidly oscillating contributions from the 
subbands of light and heavy holes are neglected, the spectrum may be both continuous and 
discrete. In the general case, near such a quasi-discrete level there exists a peak of density 
of states with small hole momenta in the direction perpendicular to the heteroboundary. The 
spectrum of holes of this kind is similar to that of light holes. The corresponding dispersion 
equation for SO holes has the form 



E g + 6 + Vc - E k 2 so + q 2 



E g + 5-E 



f^so cot /u<j cj/2 ~\~ kg 



X 



E g + 8 + V c - E k 2 so + q 



E n 



t-E K 2 so -q- 



'Kjso tctii kgod I '2 kg 



\ so E a + 8 + V c - E k 2 so + q l 



A, 



E a + 5-E 



x 



(All) 



2. Electrons 

Electrons obey the same symmetry rules as holes. Their wave functions are similar to those 
of light holes and can be written as: 

at |x| < a/2 , 

ipse — A\ cos k c x 7] + A 2 sin k c x £, 

/ k c sin k c x r] — X c q cos k c x £ \ 
ip c = —^-A\ — iq cos k c xrj + i\ c k c sin k c x^ + (A12) 

i —X c k c sin k c x £ + X c q cos k c x r\ j 

I — k c cosk c x£ — X c q sin k c xr] \ 

ih'y 

H — t^-A 2 — iX c k c cos k c x rj — iq sin k c x £ , 
\ —X c q sin k c x £ — X c k c cos k c x rj j 

(A13) 



where 



£ 2 + £{2E 9 + 26) + (Eg + 3S)E : 
£ + E„ + 25 



(A14) 



A, 



£ + E g + 28 ' 
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Here k c denotes the ^-component of the quasi-momentum of electrons in a quantum well, q 
is the longitudinal quasi-momentum of electrons. Functions with another symmetry can be 
derived by the same procedure as that used for holes. From the boundary condition follows 
that ip s and ip x must be continuous. This yields the following dispersion equation: 

fc c tan k c a/2 — ~^ K <^j ^ c cot k c a/2 + ~^ K cj = ~ <? 2 (^c ~~ ^ c ~g) ' (A15) 
where k c is the modulus of the x quasi-momentum component of electrons in the barrier region, 



Z 



£ 2 + £{2E g + 2U V + 26) + (E g + U V + 36) (E g + U v ) 



£ + E g + U v + 26 

*.= I .. 

£ + E g + U v + 26 

The spectrum splits into even and odd states if the longitudinal wave vector (g) is small or 
the expression in parentheses in the right-hand side of the equation is close to zero. The last 
condition is commonly fulfilled, since, as a rule, U v <C E g , which corresponds to semiconductors 
with about the same band structure. Note, that in the case of the discontinuous Kane parameter 
7 7^ const, continuity of "yip x and tp 8 should be usedH. 

APPENDIX B: COULOMB POTENTIAL IN THE PRESENCE OF 

HETEROBOUND ARIES 

In a quantum well the Coulomb potential of a charged particle differs from that in a ho- 
mogeneous semiconductor, because of the different dielectric constants in a quantum well and 
barrier region: 

$(r , r) = f + <£(r , r), (Bl) 

«o|r- r | 

where r is the coordinate of the particle and r is the coordinate of the point where the potential 
is observed. We consider only the case, when the particle is inside the quantum well (|xq| < a/2). 
Using the reflection method (see e.g. one can obtain: 

i.=E-f^r' (-, 1 + , 1 u 

„>i «o Vko + ^o/ \ J(x + x - (2n- l)a) 2 + p 2 J(x + x + (2n - l)a) 2 + p 2 J 




+ Y.-(^-T[ 1 + . 1 | rt|»|<./2 

^i«o\«o + «o/ \^(x -x -2na) 2 + p 2 yf(x - x + 2na) 2 + p 2 ' 

* = , 6 K0 ~ k ° + -A- £ 2n 1 + (B2) 

kqJ(x - x ) 2 + p 2 K o + K o + „>i Vk + KqJ J(x _ Xq _|_ 2na) 2 + p 2 

2C {Ko-KQ^- 1 1 

H > — , = at x > a/2. 

ko + ko^Kko + KqJ J ( x + Xo + (2n - l)a) 2 + p 2 
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Here p 2 = (y — yo) 2 + (z — z ) 2 , a is a quantum well width. These potentials are rather 
cumbersome. However, they may be simplified if dielectric constants k and k are supposed 
to be close to each other. After taking Fourier transforming along lateral coordinates (y and 
z) one obtains: 

(f)(x, q) « ( e - q \ x - x °\ + 2 /t °~ /t ° ch(g(x + x ))e~ qa ) at bl < a/2 (B3) 

2qn \ k + k J 

(f>(x, q) W — — r ( e- q{x - Xo) + K ° ~ K ° e -«(*+*o+a)\ at ^ > a/2> 

g(/t + K ) V k + /to / 

One may see that while potential itself is a continuous function along the interface, difference 
between its left and right derivatives is proportional to (kq — kq)/{kq + k ). 
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Footnotes 

1. We assume that, as is commonly the case, (V C) V v ) < E g 

2. The reason is that this term has the lowest threshold energy. By the threshold energy we 
understand the mean energy of a heavy hole taking part in an Auger transition 
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Figure Captions 



Fig. 1 Auger transition probabilities Wi(q) and w 2 (q), corresponding to the thresholdless 
and quasi-threshold matrix elements M\ and M 2 as functions of the longitudinal momentum of 
heavy holes at T = 300 K for different quantum well widths ((a) a = 50 A, (b) a = 100 A, (c) 
a = 200 A, and (d) a = 500 A). The model parameters typical of InGaAsP/InP quantum wells 
with E g ks 1 eV were used. 

Fig. 2 Threshold energy of the CHCC process as a function of the quantum well width for three 
mechanisms of Auger recombination: thresholdless (E^ h ), quasi-threshold (E 2 h ), and threshold 
(Ef h ) at T = 300 fT. The solid curve corresponds to the threshold energy of the total 

Auger coefficient (C = C\ + C 2 + C 3 ). The horizontal dashed line corresponds to the threshold 
energy Ef® for the bulk Auger process. 

Fig. 3 Auger coefficients C 1 and C 2 for the thresholdless and quasi-threshold processes as 
functions of a quantum well width at different temperatures ((a) T = 50 K, (b) T = 150 K, (c) 
T = 300 K, and (d) T = 400 K). The same parameters as in Fig. 1 were used 

Fig. 4 Coefficient C3 for the threshold Auger process as a function of the quantum well width 
at T = 150 K (a) and T = 300 K (b). 

Fig. 5 Temperature dependence of the total Auger coefficient and the partial contributions 
of the thresholdless and quasi-threshold mechanisms at different quantum well widths. The 
same parameters were used as in Fig. 1. 

Fig. 6 Comparison of the thresholdless Auger coefficient (C\) as a function of temperature at 
various Fermi energies of the holes for two different quantum well thicknesses ((a) a = 50 A and 
(b) a = 500 A). Tp denotes the Fermi energy expressed in Kelvins. The curve with Tp = —100 K 
approximately corresponds to the Boltzmann statistics. 

Fig. 7 Phonon-assisted Auger coefficients corresponding to the thresholdless and quasi- 
threshold matrix elements as functions of temperature at different quantum well widths ((a) 
a = 50 A and (b) a = 150 A). 

Fig. 8 Three-dimensional Auger coefficients C\ a 2 and (C 2 + C 3 ) a 2 as functions of a quan- 
tum well width at T = 300 K. The same parameters as in Fig. 1 were used. 
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